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Abstract 

Within the framework of the 1 /N c expansion of four-fermion interaction 
models, we analyse the next to leading 1/N C corrections to the well known 
large-iV c result Ms = 2Mq where Ms is the mass of the scalar boson 
and Mq is the constituent quark mass. The calculation is performed in 
the Extended Nambu-Jona Lasinio (ENJL) model which is suitable for 
describing low energy hadron properties. We treat the model as fully non 
renormalizable and discuss the comparison with approaches based on the 
equivalence with renormalizable Yukawa type models. We consider both 
the Gy = and the Gy ^ cases with nf = 2 flavours and study the 
dependence upon the regularization scheme. We find that pure next-to- 
leading 1/N C corrections are large and negative, while a partially resummed 
treatment can induce positive and smaller corrections. A triplet-singlet 
states' splitting is observed. 
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1 Introduction 



The physical content of four-fermion interaction models has been extensively 
analysed in the past recent years. Within the 1/N expansion approach [JTJ] for 
a general U(N) symmetric model, the equivalence under certain assumptions 
of four-fermion models with scalar four-fermion interactions and Yukawa-type 
models has been investigated in 0, |3|, while the consequences of imposing the 
so called compositness condition have been derived by [f|, [| and most recently 
by ||. In all the cases the renormalized ratio of boson and fermion masses plays 
a relevant role. In || is shown that the ratio goes to a fixed value due to the 
infrared freedom of the renormalizable Yukawa-type model and the assumption 
that the couplings are generic at the cut-off scale. In addition a common trend 
of all the analyses seems to be the fact that the large- N c value of the mass ratio 
Ms/Mq = 2 gets a large and negative next to leading 1/N C correction for a 
realistic value N c = 3. This suggests an asymptotic behaviour of the series where 
each finite order in 1/N C fails to give a good estimate of the real value of the mass 
ratio for useful values of N c . 

In this paper we address a calculation of the 1/N C next to leading correction to 
the scalar boson mass based on a treatment of a four-fermion model a la Nambu- 
Jona Lasinio which is alternative to the approaches formulated in [FJJ and || and 
is of more immediate use in the derivation of hadron properties. In this case the 
model is treated as fully non renormalizable (see M and M for reviews). 

Effective constituent quark models a la Nambu-Jona Lasinio have been found 
to be successful in reproducing the experimental values of the low-energy coupling 
constants to 0(p 4 ) in Chiral Perturbation Theory (ChPt) ||. Here the proper 
time regularization has been used. Small dependence upon the regularization 
scheme has been also verified in |]TI| [7]]. Many of the couplings between resonances 
and pseudoscalar mesons have been computed ||, |ll|] and nicely compare with 
the experiment, as well as the vector and axial- vector masses || (T2J. 

Large- N c two and three point functions have been derived in the fully fermionic 
language via the resummation of linear chains of constituent quark bubbles 
(sausage diagrams of the $ 4 theory) [Jl2| , 13 1. In all the phenomenological re- 
sults the explicit dependence upon the ultraviolet cut-off of the effective theory 
is kept treating the model as fully non renormalizable and away from the in- 
frared limit. The large- iV c calculation of the scalar two point function in the fully 
fermionic language and in the chiral limit reproduces a pole at Ms = 2Mq [12|| . 
where Mq is the constituent quark mass. With a typical value of Mq = 250 -h 350 
MeV one has Ms = 500 4- 700 MeV. The question arises if this pole has to be 
identified with a physical hadron state or it remains an artefact of the low energy 
model possibly related to the lack of confinement. One remote possibility is that 
the eventual low lying scalar resonance has a very large width. 

From the experimental point of view a signal of a narrow scalar state around 
750 MeV is reported in ||14|| , while the first clear scalar resonances are the ao(983) 
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and the isosinglet /o(975) states. Their interpretation as an ordinary qq bound 
state is dubious JL5], [L6], pTTfl . Most recently a fit of the available data |L8j indicated 
that the KK component is large for both the a (983) and /o(975) states. It is then 
clear that the identification of the physical scalar states ao and fo would probably 
require the insertion of a mixing with exotic states (glueballs, KK bound states 
etc.) inside a low energy model. This is beyond the scope of this paper. 

The present version of the ENJL model only allows for a scalar state which can 
be elementary or a qq composite state. Nonetheless we show that 1/N C corrections 
to the large- N c value of the scalar mass can produce a splitting between the octet 
and the singlet scalar states. 

In section |2| we first outline the model and make clear the main differences 
amongst the present approach and the approaches in and based on the 
equivalence of the four-fermion interaction models with a renormalizable Yukawa- 
type model. Then we clarify the correspondence between the non-bosonized 
version of the ENJL model, where only fermion degrees of freedom are present, 
and its bosonized version which only contains the auxiliary boson fields once the 
fermions have been integrated out. The appearance of overlapping divergences 
in the diagrams which give the 1/N C n.t.l. corrections to the scalar two-point 
function in the non-bosonized version can prevent from a simple and unambiguous 
calculation. We chose to compute them in the bosonized version which gives a 
reliable and fully analytical approximation of the exact result. In section |3] we 
derive the 1/N C corrections to the scalar pole mass in the Gy = case, where 
only scalar and pseudoscalar meson fields are present in the bosonized action 
and with rif = 2 flavours. Here a subsection is dedicated to the treatment of 
leading divergences in this type of theories required by chiral invariance. We also 
comment on different covariant regularization schemes. In section |4] we extend 
the model to the case Gy ^ 0, where also vector and axial-vector fields are 
present. For both cases a numerical analysis is shown and the appearance of a 
mass splitting between the scalar singlet and the non-singlet is obtained. We 
comment on numerical results and state our conclusions in section 



2 Bosonized versus non-bosonized version of the 



where C$,p and Cy s A are all the possible four-fermion lowest dimensional inter- 
actions allowed by chiral symmetry and leading in the 1/N C expansion 



ENJL model. 



The effective ENJL Lagrangian can be written as follows ||: 

£enjl = £qcd + £s,p + £v,A, 



(2.1) 
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£v,a = j^j^EMlAM^l) + (L ^ R)}. (2.2) 

c a,b 

At this level the model contains only fermion d.o.f. and is written in terms of 
three independent parameters: Gs, Gy, and the physical cut-off A of the effective 
interaction. New extra parameters can be hidden in the cut-off procedure which is 
necessary in a non renormalizable model. The couplings Gs and Gy are explicitly 
dependent upon the cut-off and we have pulled out a factor 1/N C so that they are 
0(1) in the 1/N C expansion; a, b are flavour indices and a sum over colour d.o.f 
is implicit between brackets. Cqcd ^ s the QCD Lagrangian in the presence of 
external sources and in the presence of a low energy cut-off where high frequency 
quark and gluon modes (i.e. with energy greater than A) have been integrated 
out. The problem of the connection between QCD and this Lagrangian has 
been addressed in 0. The non renormalizable (by power counting) part of the 
Lagrangian ( |2.1|) is in principle the first term of a double expansion in 1/N C , 
where N c is the number of colours, and in 1/A 2 . 

It is worth at this stage to outline the main differences of our approach with 
recent analyses of four-fermion models based on their equivalence with renormal- 
izable Yukawa- type models at least in the case Gy = 0. There are essentially 
two pictures explored alternative to the present one. A quite general RG equa- 
tions analysis of the Gross-Neveu (GN) model has been done by Zinn- Justin in 
H (see also |§ for a numerical study of a NJL model on the same lines), while 
a NJL model has been studied in ||. Here the consequences of imposing the 
compositness condition on the scalar field of the renormalizable Yukawa model as 
an additional constraint are analysed (see also refs. therein). 

In the RGE analysis in || the mass gap equation and the scalar propagator 
of the renormalizable generalized GN model in four dimensions reduce to the 
ordinary GN ones in the infrared limit a,p <C A, where a is the vacuum expec- 
tation value of the scalar field and p is the typical four-momentum. Then in this 
case the equivalence of the four-fermion model with the generalized GN model is 
strictly valid in the infrared domain within the 1/N expansion. The equivalence 
can be in principle spoiled beyond the 1/N expansion for small values of N. As 
also pointed out in || this regime can be investigated by numerical studies of 
the four-fermion model in four dimensions and compared with analytic e expan- 
sion of the renormalizable model. The type of equivalence in @ also allows for 
the presence of higher dimensions operators in the four-fermion model which are 
irrelevant in d < 4 dimensions. 

The compositness condition on the scalar field of a renormalizable Yukawa 
model treated in |J is a stronger extra constraint which guarantees the equiva- 
lence between a NJL model and a Yukawa one and which affects the RG flow of 
the renormalized couplings of the Yukawa model. It spoils the renormalizability 
of the Yukawa model in four dimensions. One underlying difference between the 
approaches in and [[]] is the fact that the compositness condition implies a 
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Figure 1: Two-point function in the non bosonized ENJL model in the large- iV c 
limit 




(a) (b) 



Figure 2: Two-point function in the non bosonized ENJL model at next-to- 
leading order in 1/N C . Diagrams (a) are self-energy insertions. Diagrams (b) 
are vertex corrections. The double lines are linear chains of constituent quark 
bubbles. 

particular value of the bare couplings at the cut-off scale while they are naturally 
assumed to be generic in the RG analysis of || . The analysis in also provides 
a prediction for the ratio of the boson and fermion masses at next-to-leading 
order in the 1/N expansion. As will be also true in our case, 1/N next-to-leading 
corrections to the mass ratio are large and negative. 

The main difference with our approach is that they keep the original four- 
fermion model at the infrared limit, which actually corresponds to the limit A — > 
oo or equivalently ff,p < A. What we do is to keep the model away from the 
infrared domain, which corresponds to being away from the limit a < A in 
the RGE analysis and in the solution of the mass gap equation. In this case 
the four- fermion model stays as non renormalizable and 1/N C corrections to the 
scalar mass can in principle be derived in the non-bosonized version where only 
constituent quarks appear. 

A full two-point function in the large-iV c limit is given by the infinite re- 
summation of linear chains of quark bubbles with the insertion of a four-quark 
interaction vertex as shown in Fig.[3]. It is of order N c (each fermion loop gives a 
factor N c ). Two and three-point functions have been derived in ||T2"| , |13|| . Analo- 
gously a full n-point function in the large- A^, limit is given by the one-constituent 
quark loop dressed by the insertion of two-point function legs attaching to the 
external sources. 
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Next to leading in 1/N C corrections are given in the diagrammatic language 
by all the possible insertions of one loop of chains of quark bubbles in the large- 
N c diagrams. They are of two types: self-energy insertions (Fig.^(a)) and vertex 
corrections (Fig.^|(b)). Being a non renormalizable model the one loop correction 
implies a new divergence and thus a new counterterm which we keep the new 
cut-off A of the loop. The exact calculation of the next to leading in 1/N C 
corrections to the scalar two-point function (and in particular to its pole mass) 
involves the one loop insertions of the type 0(a) and 0(b) in the large- N c scalar 
two-point function of Fig.|l| Defining the scalar two-point function as n(g 2 ) = 
i J d 4 x e iqx {0\TS(x)S(0)\0), where S(x) = — -j=q(x)q(x) (we omit for simplicity 
flavour indices), the large- N c expression is given by 



n(g2)( ^oo) = nm g (, 5 n(g 2 )) n = 1 _ n ^ ( ) g2) , (2.3) 

where 

S(Q 2 ) = - " (Q 2 + (2M Q ) 2 )Z S (Q 2 ) (2.4) 
9s 

is the bare fermion loop diagram in the mean-field approximation and Zs{Q 2 ) 
is the scalar wave function renormalization constant which in the proper time 
regularization is given by 

, N c r 1 ( a(l - a)Q 2 + M%\ , 

with T(0, e) = J e °° dz\e~ z and g s = Atc 2 G s /N c A 2 . 

As an example, at next to leading order in 1/N C the resummation of the self- 
energy insertion diagrams as in Fig.|(a) is given by g s U = ^+^^2^+ 
■ ■ • and one gets 

9s[ Q 2 + (2M Q p-E(QW)- 1 ]' 1 } 

with 

^ = W W? T(q ' k)2 i - «n(*») i - ss n(( 9 - *)»)■ (27) 

T(q, k) is a three-point function vertex of the type SSS, SPP, SVV, SAA, SPA 
(S=scalar, P=pseudoscalar, V=vector, A=axial- vector). 

Because of the appearance of overlapping divergences and the necessity of 
numerically evaluating contributions like fl2.7|) due to a complex momenta de- 
pendence we chose to estimate them using a reliable and fully analytical approx- 
imation within the bosonized version of the model. The correspondence with the 
non bosonized case is such that a string of quark bubbles is replaced by a meson 
line with the same quantum numbers. The n.t.l. 1/N C corrections become one 
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loop corrections in the meson theory. In the bosonized version, after integrating 
over constituent quarks, only the auxiliary boson fields remain: scalar and pseu- 
doscalar in the Gy = case and the additional vector and axial-vector fields in 
the Gy 7^ case. In what follows we refer to the bosonized version with the non 
linear realization of the chiral symmetry (i.e. the non linear representation for 
the pseudoscalar field and derivative coupling of the pseudoscalar to the other 
degrees of freedom). 

Our approximation does correspond in practice to neglecting momenta de- 
pendence in vertices and masses of ( |2.6|) . We discuss in Section [| the numerical 
relevance of the approximation. On the more formal side our approximation cor- 
responds to compute the next to leading 1/N C corrections within the bosonized 
version keeping the leading order contributions in the Heat Kernel Expansion 
(HKE) approach || to the boson vertices and masses, (see also [HJ for a review 



on HKE). Besides this it is easy to verify that the resummed HKE for a given 
interaction vertex (which is an expansion in powers of <9 2 / m q) and the large-iV c 
resummation of quark bubbles produce the SAME momenta dependence. 

For the concerns of the numerical evaluation it is useful to notice that the HKE 
behaves as a slowly convergent series (alternating signs with slowly decreasing 
coefficients) which implies that the leading term is a better estimate of the exact 
result than any truncation at a finite order outside the domain q 2 <C Mq. In what 
follows, only the lowest order in the derivative expansion will be kept for each 
vertex. This approximation allows us to simplify the calculations and preserves 
chiral invariance. 



3 The G v = case. 

In the Gy = case non renormalizable four-quark interactions in the Lagrangian 
(|2.1|) reduce to the scalar and pseudoscalar type with one coupling constant Gs- 
The bosonization introduces scalar and pseudoscalar auxiliary fields and the in- 
tegration over constituent quarks generates the effective action for the scalar and 
pseudoscalar physical mesons. The pseudoscalar sector is the ChPt Lagrangian 
of the pseudo-Goldstone bosons. The effective Lagrangian thus obtained is by 
construction globally chiral invariant (and locally chiral invariant in presence of 
external left and right handed sources) and it is non renormalizable being an 
infinite expansion in powers of derivatives acting on the meson fields. Details on 
the method can be found in 0. 

We restrict ourselves to the U(2)l X U(2)r case (we disregard the effect of 
the U(1)a anomaly which is also a next to leading effect in 1/N C ). The general 
form of the meson fields, singlets or triplets under SU(2)y, reads 

M = E ^ M (^ (a) + ^ M ol, (3.1) 
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where T^ a \a = 1, ..3 are the Pauli matrices with Tr(r a T b ) = 25 ab and Mq is 
the singlet component. In the chiral limit (m u = = 0), the effective chiral 
Lagrangian including scalar and pseudoscalar mesons at leading order 0(p 2 ) in 
the derivative expansion is given by: 

£ s > p = & < ^e>+\<d»Sd»S>- 1 -M 2 s <S 2 > + £f n r 

= -^<s*>-^<s*>+c d <s^e>+cP<s%e> 

+4 2) < S^Se > • (3.2) 



Building blocks of the Lagrangian ( |3.2|) are the scalar field S and the axial current 
of the pseudoscalar field £ M = i{^(d^ — zr M )£ — £(<9 M — il^)^}, where r M and are 
the external right-handed and left-handed sources and £ = VU = exp(--^j-) 
is the usual exponential representation with the pseudoscalar meson matrix $ 
defined as in ( [3.1|) . Both fields £^ and S transform non linearly under the chiral 
group G = U(2)l x U(2)r as O —* h{Q)Oh){$). The couplings amongst mesons 
have been derived using the HKE and with proper time regularization. Their 
expressions are listed in Appendix A. They are functions of the cut-off A of 
the fermion loop, the constituent quark mass Mq, the axial-pseudoscalar mixing 
parameter (g^ = 1 in the case Gy = 0) and the number of colours N c . 

As it is implied by the non renormalizability of the model the values of the 
parameters are a priori regularization dependent. Most suitable regularizations 
are the covariant ones: proper time, four-momentum cut-off and Pauli- Villars. 
Explicit solutions of the gap equation of the four-fermion model in the three 
cases can be found in @. A small regularization dependence of the parameters 
has been found ]7|, [10 



The next to leading in 1/N C corrections to the pole mass of the scalar two- 
point function within the bosonized version are given by the one loop corrections 
to the scalar meson propagator generated by the vertices in ( p.2|) . The diagrams 
which contribute are the ones in Fig.|3](a), the self-energy insertions and|3|(b), the 
tadpoles. The diagram |3](c) does enter the gap equation (it modifies the one point 
scalar function) and has not to be included in order to avoid double counting. 



The 1/N C corrections to the gap equation have been considered in [1C] and proven 
to be numerically relevant but still in a perturbative regime. The one boson loops 
have to be regularized and thus explicitly depend on a new cut-off A. This is the 
signal of the non renormalizability of the model. Physical inputs can be used 
to constrain its value (see section |3.2|). One loop diagrams in Fig.0 can be up 



to quartically divergent by naive power counting due to the derivative coupling 



of the pseudoscalar field. As was noticed in [20[] with a specific example of a 
phenomenological pion Lagrangian a la Weinberg, the quantization of effective 
theories like the one in (p.2[) with an arbitrary number of derivatives, can fail if 



one uses naive Feynman rules with a cutoff regularization scheme. 
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(a) 



(b) 



(c) 



Figure 3: Two-point scalar function in the bosonized ENJL model at next to 
leading order in 1/N C . Diagrams (a) are self-energy insertions. Diagrams (b) are 
tadpoles. Diagrams (c) contribute to the one-point scalar function and have to 
be included in the gap equation. 



In the next subsection we briefly show that all the leading quartic divergences 
disappear under the requirement of chiral invariance of the partition function. 
Our case is a simple extension of the example shown in |2D| 
of pions interacting with scalar fields. 



to an effective theory 



3.1 Quartic divergences versus chiral invariance 



The Lagrangian ( |3.2| ) satisfies two requisites: a) it is an infinite expansion in 
powers of derivatives acting on the fundamental fields and b) the fundamental 
fields transform non linearly under the chiral group. Expanding the field £ M 
in powers of the $ field matrix, and reducing covariant derivatives to ordinary 



ones which enter in our calculation, the Lagrangian (|3.2| ) can be written as C 



-g ($, S)d fl <& a d fl <& , where a, b are flavour indices going from (for the singlet 
case) to 3 (for the triplet case) according to the decomposition (|3.i|). The metric 
tensor g ab (§, S) is explicitly dependent on the pseudoscalar field and the scalar 
field due to the presence of interaction terms. We find 



,S6 



cab 



1 + 



4C <* r, , 2 , (1) t J2) w2 2 



+44 2) E SiSW a 5 0b + ^c?S a S b , 

J 7T i=l •* 7T 



of) 



i=l 



(3.3) 



with flavour indices a, 6=1, ...3. The metric tensor g ab defines a non linear chiral 
transformation of the pseudoscalar field $ contained in the original Lagrangian 
C = \5 ab d^ a d^^ b with a flat metric S ab . Under this transformation the full 
partition function has to be invariant (we are not concerned with anomaly in this 



S 



context). Since the chiral transformed measure is the original one multiplied by 

5 A (0)\J detg a ^, the chiral invariant partition function is defined in terms of the 

new Lagrangian C = £ + <5 4 (0) In J detg a ^. By doing an expansion in the small 
couplings to the scalar field we find 



In J detg-b = -^S, - ( -^-Y (s? + I £ StET 



i=i 



+4^ £ SiS* + y 2 (cP + cf)Sl + 0(S 3 ). (3.4) 

J 7T i=l J TV 

The new terms with <5 4 (0) = / d A kj (2-7r) 4 do exactly cancel all the leading quartic 
divergences generated by diagrams in Fig.|3| whose final expressions are listed in 
Appendix B. The first term cancels the quartic divergence of Fig||(c) (referred 
to as "top" diagram in Appendix B) for the pseudoscalar case, the second term 
cancels the one in Fig|J(a) (self-energy) pseudoscalar, and the last two terms the 
one in Fig.^(b) (tadpole) pseudoscalar in the triplet and singlet case respectively. 



3.2 Numerical analysis 

We have calculated the 1/N C corrections in two cases: 1) assuming that the 
scalar particle is a singlet and 2) assuming that the quark content of the scalar 
particle is the same as that of the p(770) vector meson. This could be the case 
of the physical ao(983) scalar resonance. Obviously our SU(2) calculation has 
to be interpreted as a first indicative approximation of the fully realistic SU (3) 
calculation. The self-energy and tadpole contributions are listed in Appendix B 
for the scalar and pseudoscalar loops and both for the singlet and triplet cases. 
In the chiral limit (m n = 0) all the pseudoscalar one loop corrections vanish. 
Denoting with M| = (2M Q ) 2 the pole mass of the scalar two-point function in 
the large- N c limit, the corrected scalar mass at next-to-leading order in 1/N C with 
a proper time regularization is the following in the singlet case: 



Ml 



Ml 



1 + T^V(-1, Ms) 
167r 



\2 
A 3 



167T 2 M| 



no, m s ) 



while for the neutral scalar triplet (the one associated with t$) we get 

1 



Ml 



Ml 



i + ?Ar(-i,i s )- 1 



\2 

A 3 



3167T 2 



2 16?r 2 M| 



m m s ) 



(3.5) 



(3.6) 



Away from the chiral limit (i.e. m n ^ 0) the additional corrections we get are as 
follows: 



Ml 2m 2 
167T 2 Jl 



Ac 2 , 

.1^(1-4 



Ml 



-1,7717, 



4(c 4 +c 4 ) M ^L 



-1,7777, 



9 



+■ 



J IT 



(3.7) 



Ml 2ml 



2c Wi _, m ' Wi m sa^)™* 



4c5 / m? 



(3i 



which also include the contribution from the wave function renormalization con- 
stant. We notice however that explicit mass terms in the pseudoscalar Lagrangian 
have not been included. The partial gamma functions T(n, Mi) of the proper time 
regularization depend upon the adimensional ratio Mf/A 2 , with i = S,P and 
where A is the new cut-off of the one-boson loop. It is worth to notice that the 
diagrams in Fig.^(c) (top) which have not been included here do not generate any 
mass splitting between the singlet and the triplet scalar component as expected 
for a contribution to the gap equation, while the self-energy diagrams give to the 
triplet component half of the contribution to the singlet one. 

We have disregarded the splitting between the singlet and the triplet compo- 
nents running in the loops. In the pseudoscalar case this is due to the U(l) axial 
anomaly, which appears in the effective Lagrangian at next-to-leading order in 
the 1/N C expansion. In the scalar case it is again a next-to-leading effect in 1/N C 
as we have shown here, although other sources can compete in this sector like 
mixing with glueballs. 

The numerical evaluation of the 1/N C corrections in Q3.5Q , ( |3.6| ), (|3.7|) and ( |3.8| ) 
needs as input the values of the large- N c parameters of the ENJL model, A and 
Mq (or alternatively A and G$) and the new one-loop cut-off A. All these quan- 
tities are regularization dependent and have to be consistently evaluated in the 
same scheme. We used the proper time regularization, while the corresponding 
expressions in the Pauli-Villars scheme can be easily obtained (previous cancella- 
tion of the spurious leading divergences due to the non invariance of the measure) 
through the substitutions listed at the end of Appendix B. For the choice of the 



numerical value of A we follow the argument developed in [ 10 1 which, although 
purely phenomenological, provides a self-consistent way of estimating the size of 
the boson loop cut-off; it proves that keeping the physical value of f n at n.t.l. 
order in 1/N C constrains the allowed range for A to be A < A max , where A max is 
of the order of the constituent quark loop cut-off A. The values for the large- A^ c 
parameters in the proper time regularization are Mq = 199 MeV and A = 667 
MeV in the Gy = case (see fit 4 of || [J). The analysis in 0j shows in addition 
a small dependence of these parameters upon the regularization scheme. In Fig.^ 
we show the squared scalar mass corrected at next to leading order in 1/N C and 
in the chiral limit (formulas Q3.5p and ( |3.6p ) in the singlet and triplet cases as a 
function of the boson loop cut-off A with fixed Mq = 199 MeV, A = 667 MeV. 



The scalar boson mass Mj in the r.h.s of ( |3.5p and ( |3.6| ) is fixed at its large- N t 
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Figure 4: The squared scalar mass at next-to- leading order in 1/N C in the case 
Gy = as a function of the ratio x = A/A for the singlet case (solid curve) and 
the triplet case (dashed curve). Here A is fixed at A = 667 MeV and Mq = 199 
MeV according to fit 4 of § 

value Mg = 4Mq. The corrections are negative both to the singlet and the triplet 
states and push the mass to zero already at A/A ~ 0.8. A triplet-singlet splitting 
is induced which grows with A but remains small. Away from the chiral limit, 
with the physical pion mass, pseudoscalar contributions are again negative but 
suppressed. The behaviour of the genuine next to leading 1/N C corrections seems 
to be in qualitative agreement with the results based on equivalence arguments 
as in [0, 0. The interesting exercise is to take into account the n.t.l. 1/N C correc- 
tions to the constituent quark mass Mq as a solution of the gap equation. This 
induces a partial resummation of the 1/N C corrections to the scalar boson mass. 
The 1/N C corrections to the gap equation have been already computed in flip" 



and they cause a positive shift of the constituent quark mass Mq as a function 
of A for fixed A and f n . In Fig.|5| we show the result of using a running value of 
Mq(A) in formulas (|3.5| ) and Q3.6Q ) which qualitatively reproduces the behaviour 



found in [ID|. The surprising result is that the partially resummed corrections 
are now positive and softer, while the splitting is not modified. 

To estimate the error which affects our zero momenta approximation we stud- 
ied the momentum dependence of each vertex entering the boson loop. All the 
couplings are weakened by the q 2 corrections and reduced in absolute value by 
about 20 -T- 30% up to — q 2 ~ A 2 . This leads to the conclusion that the q 2 re- 
summed value cannot overcome the approximated value. The same numerical 
results for the scalar mass in the Pauli-Villars regularization are obtained to a 
good approximation with the rescaling of the proper-time cut-off Apt — 2\[2Apy. 
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Singlet 
Triplet 



Figure 5: The partially resummed scalar mass squared using Mq(A) corrected 
at next-to-leading order in 1/N C in the case Gy = as a function of the ratio 
x = A/A for the singlet case (solid curve) and the triplet case (dashed curve). 
Here A is fixed at A = 667 MeV and M Q {x) = 0.199 + .0995x 2 GeV which 
reproduces a 50% of positive correction to its leading N c value at x — 1 according 
to the results in HID 



4 The G v ^ case. 

The less explored behaviour of four-fermion models is in the presence of vector 
like interactions, i.e. Gy ^ in our case. The Interaction Lagrangian of scalar 
mesons with vectors and axial-vectors at leading order in the derivative expansion 
is: 

Anf = <# } < SVuSV > +cP < S%V» > +C AP < S{^, A»} > 

+c A < SA^A" > +c [ l ] < SA^SA" > +c [ l ] < S 2 A^ > . (4.1) 

All the couplings are listed in Appendix A. Notice also the presence at 0(p) 
of the mixed term scalar-pseudoscalar-axial with coupling cap- The additional 
diagrams contributing to the scalar pole mass are again the ones in Fig.||(a), 
(b) with vector, axial, or mixed axial-pseudoscalar internal lines. All the one 
loop contributions are listed in Appendix B. In this case quartic divergences can 
be addressed to two different sources: a) for diagrams with derivative couplings 
their origin can be the breaking of chiral invariance as for the genuine pseu- 
doscalar case, b) for diagrams with non derivative couplings quartic divergences 
are a natural consequence of the the bad high energy behaviour of the massive 
vector propagator A^ u = (g^ v — k fl k u /M v )/ (k 2 — My) and they signal the non 
renormalizability of the massive vector Lagrangian. Divergences of type a) are 
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x 

Figure 6: The squared scalar mass at next-to-leading order in 1/N C in the case 
Gy 7^ as a function of the ratio x = A/A for the singlet case (solid curve) and 
the triplet case (dashed curve). The effect of using the scalar-like vector (axial) 
propagator is also shown (dot dashed curve). Here A is fixed at A = 1.16 GeV, 
g A = 0.61 and M Q = 265 MeV according to fit 1 of [§] 



cancelled following the same demonstration as in [O] where the generic field <3> is 
now replaced by a generic vector field V^. They are absent in our case. Diver- 
gences of type b) can be cured by the introduction of a spontaneous symmetry 
breaking mechanism or taking into account the compositness of the vector fields. 
Nonetheless we observe that a nearly quartic divergence could not be avoided 
in the calculation within the non bosonized version using the large N c vector 
two-point functions predicted in [P^j, where the running vector mass behaves like 
My(k) ~ In k. This is the signal of the expected bad high energy behaviour of an 
effective NJL model. In Appendix B we show the results obtained using the ordi- 
nary propagator of a massive vector field = (<7 M „ — k^/My) / {k 2 — M v ). As 
an example we also studied the results for the scalar-like propagator with softer 
renormalizable high energy behaviour A Mi , = g^ v /{k 2 — M v ). In the chiral limit 
and for the non renormalizable massive vector (axial) propagator the additional 
corrections to the scalar mass are as follows: 



AMg 1 
16tt 2 



r 2 



M 2 S 



f 2 



M 2 



2Y(-2,Ma 



1 + 



Ml 
2M\ 



T(-2,M A ) + (l - ^)r(-l,M A ) + 



T(-1,M A 



2M\. 



3- 



2M\ 



Ml 
Ml 



T(0,Mj_ 
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Figure 7: The squared scalar mass using Mq(A) corrected at next-to-leading 
order in 1/N C in the case Gy ^ as a function of the ratio x = A/ A for the 
singlet case (solid curve) and the triplet case (dashed curve). The effect of using 
the scalar-like vector (axial) propagator is also shown (dot dashed curve). Here 
A is fixed at A = 1.16 GeV g A = 0.61 and M Q (x) = 0.265 + .1325a; 2 GeV which 
reproduces a 50% of positive correction to its leading N c value at x — 1 according 



to the results in 10 



da T 



2r(-2,M A ) + 4T(-l,M A ) 



(4.2) 



AMI 






4Mf j 


' 1 


r 2 r 


16tt 2 1 


. 2 Ml . 


r 2 


Mi 


r(- 


■J IT 


Ml 



2r(-2,M A ) 



jT rf« r^o, 



+ 1 



1 + 



Ml 



1M\ 



r(-i,M A ) 



Ml 
2M\ 



2M\ 



1_M| 
,2 Mj 



aMi 



+c 



(2) My 



2 r 



MI 



A 2 

2r(-2,M y 



( 2 )Mi r 



MJ 



+ 4H 



2r(-2,M A )+4T(-l,M A ) 



r(o,M A ) 



(4.3) 



As in the = case, we studied the pure next to leading l/N c corrected 
scalar mass and the partially resummed one using for the vector and axial masses 
My = 0.8 GeV and M A = 1 GeV. The first result is shown in Fig|| The large- 
N c values of the parameters in the Gy ^ case with proper time regularization 



are M Q = 265 MeV, A = 1.16 GeV and g A = 0.61 (see fit 1 



in 



A 
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few comments are in order. The singlet-triplet splitting is enhanced respect to 
the Gy = case. The singlet mass still receives negative corrections. The 
anomalous enhancement of the triplet mass is sensitively dependent on the form 
of the propagator. It is consequence of the presence of positive contributions 
in the vector sector which actually dominate in the case of the ordinary vector 
propagator form and that are zero in the singlet case. The partially resummed 
behaviour is shown in Fig.|7|. Corrections to the singlet state are softened but still 
negative. The same anomalous enhancement of the triplet mass is observed. In 
both cases the singlet-triplet splitting is enhanced respect to the Gy = case. 
Again, on the base of the study of the q 2 dependence of the vector (axial) vertices 
we expect that the inclusion of the full q 2 dependence will soften the corrections. 
Within the present approximation the largeness of the axial and vector corrections 
prevents from a fully reliable estimate in the region A/A ~ 1. 

5 Conclusions 

We studied the next-to- leading in 1/N C corrections to the pole mass of the scalar 
two-point function within the bosonized version of the Extended NJL model and 
away from the infrared domain. In this context the model is treated as fully 
non renormalizable and a new cut-off parameter have to be introduced for the 
one boson loop. Within a reliable zero momenta approximation, which is the 
leading order of the Heat Kernel Expansion, we have analytically derived the 
next-to-leading 1/N C corrections to the scalar mass in both the Gy = and 
Gy 7^ (vector and axial fields present) cases and studied their regularization 
scheme dependence. The main results are that genuine next to leading 1/N C 
corrections to the singlet state are negative and relatively large, while a partially 
resummed estimate induces positive and softened corrections in the Gy = case. 
Remarkably the corrections to the large- N c degenerate mass for the triplet and 
singlet states induce a splitting which mimics the physical one (octet heavier than 
the singlet). The splitting effect is enhanced in the Gy ^ case. 

The largeness of the negative pure next to leading 1/N C corrections derived 
in this framework qualitatively agrees with the results for the scalar over fermion 
mass ratio derived in the IR limit where the equivalence with renormalizable 
Yukawa- type models is valid by the use of the compositness condition ||. This 
suggests an asymptotic behaviour of the 1/N C expansion of the mass ratios in 
four- fermion models both at the IR limit and away from the IR limit, where a 
truncation at any finite order fails to be a good estimate of the real value for 
useful values of N c . 
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A The couplings of the bosonized Lagrangian 



The scalar-pseudoscalar couplings are: 



A3 
3! 

A4 
4! 



c 4 



J2) 
C 4 



N r .Mr 



167T 2 ZT 



r(o,e)--r(i,e) 



Nr. 



16tt 2 Zl 



r(0,e)-4r(l,e) + -r(2,e) 



-Mr 



lQn 2 "^y/Zs 
1 N c g\ 



2 16tt 2 Z s 
1 N c g\ 



r(o,e)-r(i,e) 



20. 



2 16tt 2 Z s 

The scalar-vector and scalar-axial couplings are 



r(0,e)--r(l,e) + -r(2,e) 
r(0,e)-^r(l,e) + ^r(2,e) 



1 N r . 



M, 



ca 



C A 



C A 



cap 



Q 



2 16tt 2 JZ~ s Zv 



4r(0,e)-4r(l,e) 



1 Nr 



1 



2 16tt 2 Z S Z V [ 

1 _Nc L 

2 16tt 2 Z s Z x 



r(0,e)-^r(l,e) + ^r(2,e) 
on c 

r(o,e)-^r(i,e) + |r(2,e) 



1 N c g A M Q 

2 16tt 2 ^[Z~s^fZ^ 



-4r(0,e)+4r(l,e) 



-V 



-c (2) 

Cy 



1 N r 



(A.l) 



2i6. 2 ^r r(0 ' e) + i r(M) J- (A - 2) 

All the couplings have been derived within the Heat Kernel Expansion with 
proper time regularization. Zy and Za are the wave function renormalization 
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constants of the vector and axial- vector fields Zy = Z A = N c /A8n 2 T(0, e). 
The partial gamma functions T(n — 2,e), with e = Mq/A 2 , are defined as 
T(n — 2,e) = / e °° dzl/ze~ z z n ~ 2 . T(— 2, e) contains a quart ic divergence, T(— l,e) 
a quadratic one, T(0, e) is logarithmically divergent, while T(n, e) with n > are 
finite. 



B The one loop contributions 

Here the contributions of the three classes of diagrams of Fig|J are listed for the 
scalar (S), pseudoscalar (P), vector (V), axial (A) or mixed axial-pseudoscalar 
(A-P) bosons running in the loop. They are the self-energy diagrams of Fig.|3](a), 
the tadpole diagrams of Fig.|3](b) and the top diagrams of Fig.|3](c). In the case 
of axial and vector loops we give the result for the form of the propagator = 
(g^ - k^/M v )/(k 2 - M 2 ) and the softer one = 9lMU /(k 2 - My). The 
self-energy contributions are written in the form A + Bq 2 with q 2 Minkowskian 
and where B gives the wave function renormalization constant which enters the 
correction to the scalar mass. 



Self-energy diagrams 
Singlet propagator 



S 
P 



A (A 



a (Ay 

A-P(A^ U ) 



AC 2 /2^r 



16vr 2 V/ 2 



-2m*r(-l, m n ) + m^r(0, m n ) 



-I 2 ( ^^(-1, m n ) - m 2 r(0, m n ) 



7? r 



4c^ 
16vr 2 

2M\ 

7,1 



2r(-2, m a ) + r(-i, m a ) + 3r(o, m a ) 



r{-l,M A )-T{0,M A/ 



451 
16vr 2 



4r(0,M A ) 



167T 2 / 2 



I 3 TTt 

-M 2 r(-2,M A ) + -^r(-2,m.) 



K + M|) 2 \ 
AM\ 



da T O 



[I - a)ml + aM 2 A 
A 2 - 
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1 



ml + Ml 



2M\ 



a)m\ + aM\ 



A-P(A') = i 



8c 2 AP 2 



167f 2 f 2 



£ da r(o,± 

M\Y{-1, M A ) - ml J* dav(o, 



a)m\ + aM\ 



A 2 



(B.l) 



Notice that no vertex SVV is allowed. 

Triplet propagator 

The only possible self-energy diagrams for the triplet propagator contain two 
different internal lines, one singlet and one triplet. The contribution is half the 
contribution to the singlet propagator displayed above. 

Tadpole diagrams 
Singlet propagator 



S 
P 



A, 



167T 2 

4(c?> + cf 



Mfr(-1,M S ), 



167T 2 fl 



m^r(-l,m 7r ) 



A(V)(\„) = 



a(v)(a; u ) = 



4(c 



(i) 



S 2 ) 



16tt 2 



4(c' 



(i) 



„( 2 ) 



A( V) ^A(V)j 

16tt 2 



2T(-2,M A{V) )+AT(-1,M A{V) ) 



4F(-1,M A{V) ) 



(B.2) 



Notice that in the singlet case there are no contributions from vector vertices of 
order O(p ) like < S 2 V^V^ >, because Cy ] + Cy ] = 0. 



Triplet propagator 



S 
P 

MY) (A^) 



.2 A 4 



3 16tt 2 



M|r( 



-1,M S 



4c 



(i) 



167T 2 / 2 



m^r(-l,m 7r ) 



4c 



— «• 



(2) 



Win 2 ^) 



2r(-2,M A(v) )+4r(-l,M A(v) ) 
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Top diagrams 



—i 



4c (2) 



4F(-1,M A{V) ) 



(B.3) 



5 
P 



A f A 



A (AM 



x 2 



16tt 2 V ' A; 
2A 3 c d 2 1 4 



—2 



1 16tt 2 Mj 
. 2A 3 c A Ml 
* 16tt 2 Mj 



2r(-2,M A ) + 4r(-i,M A ; 
4r(-i,M A ) 



(B.4) 



Contributions are the same for singlet and triplet propagator. 



Pauli Villars Regular izat ion 

For the comparison with the proper time regularization contributions the 
following substitutions can be performed: 



m 2 r(-l, m) -> m 2 [(l + 2x) ln(l + 2x) - 2(1 + x) ln(l + x) 
T(0, m) -> 2 ln(l + x) - ln(l + 2x), 



(B.5) 



where x = A 2 /m 2 . This corresponds to the usual Pauli Villars procedure in a 
scalar theory where two additional fields with masses Mi = m + A and M 2 = 
m + 2A and coefficients Ci = —2 and 6*2 = 1 are sufficient to make the theory 
finite. In the case of a non linearly realized symmetry (as in this case for the 
pseudoscalar sector) quartic divergences due to the non invariance of the measure 
have to be treated before. 
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